Interaction-induced impeding of decoherence and anomalous diffusion 
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We study how the interplay of dissipation and interactions affects the dynamics of a bosonic 
many-body quantum system. In the presence of both dissipation and strongly repulsive interactions, 
observables such as the coherence and the density fluctuations display three dynamical regimes: an 
initial exponential variation followed by a power-law regime, and finally a slow exponential con- 
vergence to their asymptotic values. These very long-time scales arise as dissipation forces the 
population of states disfavored by interactions. The long-time, strong coupling dynamics are under- 
stood by performing a mapping onto a classical diffusion process displaying non-Brownian behavior. 
While both dissipation and strong interactions tend to suppress coherence when acting separately, 
we find that strong interaction impedes the decoherence process generated by the dissipation. 
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Understanding the influence of the environment on the 
dynamics of physical systems is of paramount importance 
in the development of quantum-based technologies. For 
systems as diverse as those encountered in solid-state 
physics [l|-|4, atomic and molecular physics [5[ or quan- 
tum optics p, lZ| 1 coupling to an environment often re- 
sults in the loss of coherence. Devising methods to con- 
trol and suppress decoherence is therefore at the heart 
of various present and future applications such as high- 
precision clocks, measuring devices [8|, |9j, and quantum 
computers [lO| . A clear understanding of the effects of an 
environment on correlated quantum many-body systems 
is presently lacking. For example, the consequences re- 
sulting from the interplay of interactions and dissipation 
are not well known. Even though each process individ- 
ually tends to suppress coherence, it is not understood 
if they cooperate or destructively compete when acting 
simultaneously. An ideal tcstbed to answer such ques- 
tions is provided by cold atomic gases as, in these sys- 
tems, both the interactions and dissipative processes are 
highly controllable. 

Exciting phenomena have already been proposed to 
occur when cold gases are coupled to an environment: a 
Zeno-like behavior due to local atom losses in systems 
of interacting atoms confined to double wells |lll4l3 or 
optical lattices [14| ; improved stability against weak dis- 
sipative effects for bosonic atoms in both the weakly in- 
teracting 15| and Mott-insulating regimes [16| ; a dynam- 
ical phase transition between a condensed and a thermal 
steady state [iTJ for a coherence-enhancing dissipative 
process; and the possibility to eng ineer dark states with 
highly desirable properties |17h22| . 

In this work, we investigate the influence of a Marko- 
vian (memory-less) environment, e.g. a light fleld or a 
thermal cloud, on strongly interacting bosonic atoms in a 



double well potential. We find that dissipation and inter- 
actions destructively compete, resulting in two surprising 
phenomena: (i) contrary to naive intuition, decoherence 
is slowed down in the presence of strong interactions, (ii) 
at intermediate times, physical quantities such as the co- 
herence and density fluctuations are found to display a 
slow power-law time-dependence characterized by a time- 
scale t* which is considerably enhanced by interactions. 

The algebraic regime can be characterized in detail 
in the limit of a large number of atoms N by mapping 
the evolution onto a classical diffusion in the configu- 
ration space of all Fock states. For A^ bosons confined 
to two wells, this space is described by a single coordi- 
nate X. The diffusion rate is found to be strongly de- 
pendent on the configuration: it is large (small) for con- 
figurations which have a low (high) cost in interaction 
energy. As shown analytically, this translates into a non- 
Brownian diffusion [23| with \/{x'^) ^ t^'^ and a non- 
Gaussian probability density p{x,t) cx r^^'^e"^ '^'^. 
Here r = t/t*, with the long time-scale t* given by 
7t* = 2N'^{U/J)^ (7 is the Markovian dissipation rate, 
U the interaction strength and J the tunneling coefficient 
between the two wells). 

Considering N bosonic atoms trapped in a sufficiently 
deep optical lattice with L sites, the evolution of the sys- 
tem is described by the master equation 



ihdtp = [H,p] +ihV (p) . 



(1) 



The first term [H,p\ describes the unitary evolution 
of the density matrix p governed by the Hamiltonian 



H = -JEti 



h\hi+i + h.c)+^Y.ti f^i {f^i - !)• The 



operators h\ and bi arc bosonic creation and annihi- 
lation operators on site / and fii = bjbi counts the 
number of atoms. At large U/J the atoms tend to 



be more localized which typically results in a reduc- 
tion of the coherence between neighboring sites C = 
Y^i^fii+i + h\j^ihi). The dissipative process is modeled 

by I?(/5) = 7Z];=i {nipni - \nfp~- \pnj), where 7 is 
the coupling to the environment and the quantum jump 
operators correspond to local density operators fii. This 
form of dissipation has been identified as one of the most 
important heating processes in the presence of a red- 
detuned optical lattice potential [16|, |2J| . In the absence 
of interactions [U = 0), it causes a rapid exponential loss 
of coherence C{t) = C(0)e-T*. 

It is useful to recast the quantum master equation ([l]) 
into an eigenvalue equation Mp = — Ap, with p^ and 
— Aq being respectively eigenvectors and eigenvalues of 
the matrix M. The positive real parts of Aq. are the in- 
verse relaxation times, since Pa{i) — s~^°'*Pa{^)- Hence, 
the longest living states are associated with the smallest 
Re(AQ). For J = 0, all combinations of diagonal Fock 
states |{n-i})({'T-/}| are steady-states (Aq = 0). However, 
for J 7^ (and L finite), the unique steady-state is the 
completely mixed state p^ ex X]{n } l{'"-'})({"'}l i^ which 
all Fock states have equal weight (corresponding to max- 
imal Von- Neumann entropy). We will show that p^. is 
reached in a highly non-trivial way. 

To pinpoint this, we concentrate on the two site Bose- 
Hubbard model {L ~ 2). Various aspects of the in- 
fluence of dissipation have been studied in this set-up, 
focusing mainly on weak interactions 25l429|. For our 
study the atoms are prepared in the symmetric ground 
state of the Hamiltonian and the dissipation is turned 
on at time i = [30|. We monitor different experi- 
mentally measurable quantities, such as the coherence 
between the two wells C, the local density fluctuations 
K = {n\) — (ni)^, and the probability P^ to measure a 
balanced configuration with N/2 particles in each well. 
At low interaction strength U or small particle number 
TV, exponential evolutions dominate. Their time-scales 
depend crucially on the parameter regime considered. In 
contrast, for sufficiently large N and UN ^ ^7, three 
distinct regimes can be identified: (i) a fast exponential 
variation at short times (7^ ^1), (ii) a drastic slowing 
down associated with a power-law regime at intermedi- 
ate times (1/7 ^ t ^ t*), and (iii) a slow exponential 
approach to the asymptotic value at long times {t >t*). 
These regimes are exemplified in Fig. [1] for the rise of 
density fiuctuations and the decay of coherence obtained 
by solving Eq. ([1]) numerically. 

A first insight into the origin of regimes (ii) and (iii) can 
be gained by analyzing the eigenvalue spectrum (inset of 
Fig. [ija)). The final exponential regime is due to the 
fact that the lowest non- vanishing real part of an eigen- 
value A (X jJ^/U^N^ is isolated. The algebraic regime 
arises from the complex interplay of many slowly decay- 
ing states Pq, (having non-zero overlap with the initial 
state) corresponding to the band of eigenvalues in the 
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FIG. 1. (color online) (a-b) Density fluctuations (a) and 
rescaled-coherence (b) versus rescaled time r. A data collapse 
is observed for large interaction strength U and atom num- 
ber N, independently of the other system parameters. The 
analytical results are taken from Eq.((6]) and ((Tjl. The inset 
of (a) shows the real parts of the smallest Xa versus N and a 
fit oc 1/N^. (c-d) Coherence decay at short (c) and long (d) 
times (A'^ = 60, U / J — 20, and hy/J = 1) with exponential 
fit (continuous green line). L = 2. 



inset of Fig.^a). 

To understand quantitatively the different dynamical 
behaviors, we reduce the system to a simpler model for 
t ^ 1/7 and strong interaction U ^ J, h'y. This in 
turn can be mapped using the large N limit onto a clas- 
sical diffusion problem in configuration space. We rep- 
resent the density matrix in the basis of Fock states as 
p = ^ p„^„i|7i)(77i| where n,m = 0,1,..., iV labels the 
number of atoms in the left well. In this regime, it is 
justified to use adiabatic elimination and map the full 
evolution onto an equation for the diagonal elements p„.„ 
only. This is done by integrating dtPn,n+i (derived from 
Eq.©) by noticing that the diagonal terms of the density 
matrix arc slowly varying compared to the off-diagonal 
ones (following :31j and for even A^). This gives 



iW^JrIi 



ih'-f 



where Jr = 2j7JV^ 



Wn 



{N - 2n - 1)U 



_ (n+l){N-n) 



Apn{t) (2) 



nHp 



'n+l.n+l 



Pn 



"+1 " (Ti-Af/2+l/2)2 ' ^l^*^ ^Pn — 

). Using ^ into dtPn.n, one obtains 



drPn.n = 2 (W,j+iAp„ - W„Ap„_i) . 



(3) 



Remarkably, the interaction, dissipative coupling, tunnel- 
ing amplitude and particle number only enter the equa- 
tions for the diagonal elements p„.„ via the dimension- 
less time T = t/t* . The emerging time-scale t* = "^^ Ji 
controls the long-time behavior of the system, and be- 
comes very large for strong coupling U and large N |32l | . 
Additionally, this dependence explains the scaling of the 
lowest eigenvalues. 

We now perform the large N limit calculation in which 
the discrete master equation ^ is mapped onto a classi- 
cal diffusion equation in the continuum. We map the 
configuration space onto a coordinate space by x = 
n/N — 1/2 e [—1/2,1/2] which becomes continuous in 
the large N limit. The boundaries {x ~ ±1/2) corre- 
spond to strongly imbalanced occupation of the double 
well, whereas the center {x = 0) corresponds to the bal- 
anced configuration. The diagonal elements of the den- 
sity matrix are related to a probability density p(x,r) by 

1 /9 

Npn.niT) = p{x,t), with normalization J_.,„p{x)dx ~ 1 
(insuring tr/5 = 1). We thus obtain the diffusion equation 



dx [D{x)dxPix,T)] = dr p{x,t). 



(4) 



The diffusion function D{x) — -^ — 1 is strongly depen- 
dent on the variable x: it diverges at the center (x = 0) 
and vanishes at the boundaries (x « ±1/2). Physically 
this slow diffusion at the boundaries corresponds to the 
slow population of the energetically costly imbalanced 
configurations. 




FIG. 2. (color online) Diagonal terms of the density matrix 
p-n,n versus rescaled time r in log-log form: the exact solu- 
tions of Eq. IT]) for n = 41(top)...80(bottom) are represented 
by blue continuous lines; the element n = 40 (thick red line) 
is proportional to the probability of finding a balanced config- 
uration Pi,; the corresponding difltusion density p(x = n/N,T) 
Eq. ([5]) up to n/N = 60 are shown in dashed green lines. Inset: 
3D plot of the evohition. Parameters: U / J = 20, h'^ / J = 1, 
iV = 80 and L = 2. 

Within this mapping, the initial ground state corre- 
sponds to p[x,T = 0) and is peaked at x = 0. The 
asymptotic long time limit is the uniform distribution 
p{x,T = oo) = 1 representing the totally mixed state 



PS- To gain more insight into the actual diffusion pro- 
cess, we insert the scaling ansatz p{x,t) = —ufiO with 
£_ = x/t^ into (HI). We find that a scaling solution exists 
provided z/ = 1/4 and r ^ 1. The function /(^) satisfies 
^/" + (C'' - 2)/' + CV = 0, whose solution can be found 
in closed form: /(^) ex exp(— ^*/4). Therefore, the diffu- 
sion process at short rescaled time (r ^ 1, i.e. t <C i*) is 
non-Brownian and described by 



p{x,t) 



V2 



1 



r(l/4) ri/4 



exp I 



-a;V4r) 



(5) 



where r(l/4) in the normalization denotes the gamma 
function. Compared to normal diffusion given 

by p{x,t) = l/\/4:TT(TT exp(— x^/4(jt) for a space- 
independent diffusion constant cr, the divergence in D{x) 
for a: « leads to a highly accelerated initial diffusion. 

The analytical expression ([5|) reproduces accurately 
the exact evolution of the diagonal terms of the density 
matrix computed from ([1]) as shown in Fig. [2] The very 
fast initial broadening of the peak is followed by a much 
slower occupation of the boundary (imbalanced) regions. 
In particular, we observe a clear power-law decay (with 
exponent 1/4) in the exact evolution of the density ma- 
trix elements, in very good agreement with the analytical 
form ([5]) even for a relatively low number of atoms. For 
times T = t/t* ^0.1 (time-scales which can actually be 
very long in an experiment), the lowest eigenvalue of the 
discrete spectrum of Q due to the finite atom number 
controls the eventual exponential convergence of p{x,t) 
to the uniform distribution. 

The diffusion in configuration space has important con- 
sequences for experimentally measurable quantities. The 
probability of measuring a balanced configuration Pb{t) 
is directly proportional to the central density p{x = 0,t). 
Therefore, its detection would reveal the passage from the 
initially exponential decay to the algebraic one in t~^'^ 
(sec the red thick line in Fig. [2]). The density fluctuations 
K show an increase which in the power-law region is well 
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described, using k/N'^ = /_-^ ,2(^^ + x)p{x,T)dx, by 



K _ 2r(3/4) 



r(i/4) 



V2r(3/4) J 
r(l/4) NU 



(6) 



As shown in Fig.[T](a), this equation describes accurately 
the exact density fluctuations at intermediate times with- 
out any adjustable parameter. At longer times, the 
density fluctuations tends towards its asymptotic value 
■^ ^ 12 + QN- For the coherence, from the continuum 
limit of ([2]), we can derive a scaling behavior C{t) = 
J/UC{t/t*) where C(r) « /^^(J^ ^^-^i^dxp{x,T)dx. This 
scaling is verified over a large range of different parame- 
ters t/ > fi.7 , J as shown in Fig. [Ub). In the power- law 
regime, where r is small, this reads 



C{t) _ r(3/4) J 1 



r(3/4) 



1 



N 



2r(l/4) UN ^ V2r(l/4) ^/7t 



(7) 



with again excellent agreement with the exact results 
over a large range of intermediate times without any ad- 
justable parameter (Fig. [TJb) and Fig. [J). We note that 
in this regime the coherence only depends on the original 
dissipative coupling 7 and not on J or ?7. 
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FIG. 3. (color online) Coherence divided by number of atoms 
versus time for different interactions and for N — 60, h"//J = 
1 and L — 2. Inset: zoom on the short times. 

In Fig. [3] we compare the decay of the coherence for 
stronger and weaker interaction strengths. A weaker 
interaction strength typically leads to a larger value of 
the coherence in the initial state and at small times (in- 
set). After the initial exponential decay, the regime of 
algebraic decay ex l/-y/t is manifest for large interaction 
strength UN ^ Hj, J. This regime shrinks when lower- 
ing the interaction, and is followed by the eventual expo- 
nential decay at long times. Surprisingly, the coherence 
for larger interactions exceeds that for lower interactions 
and survives for much longer times. Hence the presence 
of interactions impedes the dissipation-induced decoher- 
ence. The physical reason behind this finding is that it is 
difficult to populate Fock space configurations which arc 
energetically disfavored by strong interactions. We stress 
that the slowing down of decoherence, unlike the use of 
particular squeezed, "decoherence- free" , states |33[ , gen- 
erally holds independently of the initial condition. 

To summarize, we found that the time-evolution of 
the coherence and density fluctuations of strongly inter- 
acting bosonic systems subject to dissipative effects can 
present three consecutive regimes: an initial fast expo- 
nential one, and an intermediate power-law regime asso- 
ciated with anomalous diffusion which is eventually fol- 
lowed by a final exponential regime. The latter sets in 
after a time ^ t* which becomes very large for strong 
interactions and large number of particles, and the decay 
rate is reduced accordingly. Simulations on three-site 
systems and mean-field calculations for an infinite lattice 
confirm the existence of the slowing down of decoherence 
and the presence of a power law regime. We understand 
that a power-law dynamics, atypical for the considered 
Markovian environment, also occurs for different kinds of 



To probe these phenomena experimentally, one could 
use a single double well potential or an array of decou- 
pled double wells generated by an optical super-lattice 
potential. Observing the different regimes proposed here 
requires strong interactions while maintaining the sin- 
gle band approximation during the whole experiment. 
To meet this requirement, we suggest the use of a light 
species, for example Lithium, whose interaction can be 
tuned by a Feshbach resonance [35|, l36| . The light mass 
enables the use of a deep optical lattice potential (a band 
gap of about 250 kHz is realistically achievable) while 
keeping reasonable time-scales. An alternative realiza- 
tion would be to trap two bosonic species in one well and 
to tune their interaction with a Feshbach resonance [37 1 . 
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